Localized Modes in Open One-Dimensional Dissipative Random Systems 
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We consider, both theoretically and experimentally, the excitation and detection of the localized 
quasi-modes (resonances) in an open dissipative ID random system. We show that even though 
the amplitude of transmission drops dramatically so that it cannot be observed in the presence of 
small losses, resonances are still clearly exhibited in reflection. Surprisingly, small losses essentially 
improve conditions for the detection of resonances in reflection as compared with the lossless case. 
An algorithm is proposed and tested to retrieve sample parameters and resonances characteristics 
inside the random system exclusively from reflection measurements. 
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Random stratified media are found in numerous geo- 
logical and biological settings as well as in fabricated ma- 
terials. Wave interactions in such systems determine, for 
example, reflection and transmission from multilayer di- 
electric stacks used as optical reflectors, filters and lasers, 
and propagation of seismic waves in the earth's crust, mi- 
crowave radiation in sand layers, and sonic waves in the 
oceans. In addition, though considerable effort has been 
expended to develop highly periodic structures, devia- 
tions from periodicity can significantly modify the char- 
acteristics of optical and microwave tunable photonic 
crystals. At the same time, it may be possible to utilize 
highly disordered samples, for many applications. For 
example, tunable switches or narrow line laser sources 
[H can be created in randomly stacked systems or in dis- 
ordered fibers. Since the location, width, and intensity 
of resonance are random, these can be found only by di- 
rect measurement of the field inside the sample, which 
is generally not feasible. Although strong localization of 
waves and resonances in one-dimensional random media 
has been extensively studied theoretically @, H, 0] , most 
of the analytical results were obtained for lossless sys- 
tems and for values averaged over ensembles of random 
realizations. But results for an ensemble do not provide 
the intensity spectrum within a particular sample, which 
is often essential. 

In this letter, we consider samples, in which dissipation 
is weak in the sence that the condition 



—h oc = Yli oc < TL < 1, 



(1) 



is satisfied. Here, T t and T, are, respectively, the tem- 
poral and spatial decrements of the wave energy due to 
loss, v g is the local group velocity of the wave, L is the 
sample length, and h oc is the localization length. 

The nature of wave propagation in absorbing random 
systems is illustrated by measurements in a single-mode 



microwave rectangular waveguide filled with random el- 
ements. Samples are composed of random mixtures of 
31 randomly oriented elements which are 8 mm in length 
and are comprised of a solid ceramic first half and a sec- 
ond half which is milled to create an air space between 
two thin ceramic walls of thickness 0.8 mm, five 4-mm- 
long, high dielectric material slabs, identical to the first 
half of the binary element, and five 4 mm-long, low di- 
electric Styrofoam slabs. Only with the introduction of 
the 4 mm elements are states introduced deep within 
the band gap. The samples are prepared by choosing a 
random order of the elements and a random orientation 
of the binary element (ceramic part facing eather input 
or output) in the waveguide. Comparison of measure- 
ments of the frequencies of the first mode at the two 
band edges of the band gap in a periodic structure of bi- 
nary elements with a one dimensional scattering matrix 
simulation which includes waveguide dispersion gives nl 
~1.67 and n2 ~1.08 for the refractive indices of the first 
and second portions of these elements, respectively. The 
field inside the waveguide is detected along a 2 mm wide 
slot along the length of the top of the waveguide. The 
field is weakly coupled to an antenna so that it is not 
affected by the measurement. The core of the coaxial 
adapter without a protruding antenna picks up the field 
inside an enclosure coupled to the waveguide through a 
2 mm hole which is pressed against the slot. For each 
random configuration, the detector is translated in 1 mm 
steps along the slot. The in and out of phase components 
of the field were collected at each position with use of an 
HP8720C vector network analyzer. Spectra over 6 GHz 
were collected in 3.5 MHz steps. The cutoff frequency of 
the empty waveguide is 13.5 GHz. The frequency range 
covered is slightly larger than the stop band of the peri- 
odic structure. 

Figure 1 shows the intensity, I(x, v), inside two random 
configurations normalized to the intensity of the incident 
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FIG. 1: Intensity VS frequency and position inside the sample 
for two different samples. 



wave as a function of coordinate and frequency. Distinct 
localized states with high intensity are seen in Fig. la. In 
addition, multiple peaks in space such as those shown in 
Fig. lb are found in other configurations. This field dis- 
tribution has been discussed theoretically by Lifshits and 
Kirpichenkov and by Pendry || . Such a distribution has 
been termed a "necklace state" , though they correspond 
to several spectrally overlapping states rather than to a 
single state @. These states are associated with non- 
Lorentzian spectra, which have been observed in optical 
transmission through random layered media [a, lZ( . 

The reflectance spectrum, R(v), for the distributions 
Fig. la is shown in Fig. 2a. The instrumental error of 
the measurements of the reflection coefficient was about 
5%. The transmittance (intensity at the output), T(v) = 
I(L,i>), within a spectral range which encompasses only 
exponentially peaks modes on resonance and evanescent 
waves off resonance is below the experimental noise level 
(Fig. 2b). Sharp minima are observed in R(v) that are 
not accompanied by an essential increase of T(v) (note 
the difference in the scales in Figs. 2a and 2b). This is in 
contrast to the lossless case where R and T are connected 
by energy conservation, R + T = 1. Dips in R{v) arise 
because the energy lost at each point within the sample 
is proportional to the product of the intensity of the field 
at that point and spatial absorption coefficient. Energy 
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FIG. 2: Spectra of (a) reflection, (b) transmission, and (c) 
normalized energy loss, Ii OS s{v), for Fig. la. The reflection 
spectrum of the sample with the distribution shown in Fig. lb 
is shown in (d). 



conservation requires, 



1 — T — R = I loss = T / -2-1 (x, v) e{x)dx, 

JO VgQ 



(2) 



where Ii oss is the rate of energy dissipated in the sample, 
e(x) is the dielectric constant, and v g o is the group ve- 
locity in the empty waveguide. For the following estima- 
tions, permittivity e(x) is replaced by its mean value e, 



whereas v g o ~ c and v g ~ c/s/e. Then, the rate of energy 
dissipation can be written as Ii oss ~ T\fi I (x, v) dx. 

Off resonance, the ensemble average intensity falls ex- 
ponentially with decay length i; oc . The energy loss can 
therefore be estimated as Iioss ^ r^/ocV^ 1* Since the 
off-resonance transmittance T ~ T t y P = exp {—L/li oc ) <C 
1, the reflectivity R is close to unity, as it would be 
in the lossless case, but it saturates at R ~ Rt yp = 
1 — r^ocv 7 ? - T typ . 

On resonance, loss can be of order unity due to the 
high intensity in localized states, even when Eq. |1| holds 
(Fig. 2c). This suppresses both the reflected and trans- 
mitted fluxes. In Fig. 2, weak loss is seen to suppress the 
transmitted wave on resonance below experimental noise, 
while the suppression of the reflected wave is easily ob- 
served. The "necklace" shown in Fig. lb is also easily 
detected in reflection as a group of closely spaced peaks 
with the number of peaks equal to that of the coupled 
states underlying the "necklace" (Fig. 2d). 

To interpret the experimental data presented above, 
we extended the approach developed in [§] for modelling 
lossless random systems by incorporating small absorp- 
tion in the model. The approach is based on associating 
each localized state at a frequency v with a structure 
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comprised of an essentially transparent segment ( "well" ) 
of small length I <C L, surrounded by essentially non- 
transparent segments ( "barriers" ) with small transmit- 
tance T 12 *C 1 characteristic of "typical" non-resonant 
configurations. At the frequency v the configuration 
forms a high-Q resonator and wave propagation through 
it can be treated as a particular case of the general prob- 
lem of the transmission through an open resonant system, 
regular or random, whether it is quantum-mechanical po- 
tential well, optical or microwave resonator or ID ran- 
dom medium. The distinguishing feature of the random 
medium is that there are no regular walls, and the sepa- 
rate transmission coefficients of the isolated left and right 
segments of the sample, T\ and T 2 , respectively, are ex- 
ponentially small as a result of Anderson localization. 

A necklace state exists in a random configuration when 
at some frequency there are several effective cavities sep- 
arated by typical non-transparent segments . Generaliza- 
tion of the theory to this case amounts to the solution of 
a bit more complicated problem of the resonant propaga- 
tion through a set of spatially separated potential wells 
(coupled resonators) whose number is equal to the num- 
ber of dips in the measured frequency dependence of the 
reflection coefficient (Fig. lb). It is beyond the scope of 
this letter and will be considered in a future publication. 

Calculations of the resonant transmittance and re- 
flectance of an open resonator with the absorption rate 
r, and of the intensity inside the effective cavity yield 
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riVir res /T 2 , 



^res — 



(3) 



(4) 



(5) 



The half-width of the resonant peaks in R (y) is given by 



Sv 



■T1+T2 



(6) 



In the case of Anderson localization in a disordered 
sample, coefficients T1.2 are approximately given by 
Ti : 2 — exp [— (L/2 ± d) /Z/ oc ], where d is the coordinate 
of the resonance relative to the center of the sample, 
so that T X T 2 ~ T typ . Eqs. ©-© satisfy the conser- 
vation law Eq. |2]) if Ross — I r esFl\/I. The reflection 

coefficient, R res , and the rescaled intensity at the res- 
1/2 

onance, I* es — T t 'I res , depend on the position of the 
cavity and on the absorption through universal functions 
of two dimcnsionlcss parameters a = exp (— d/h oc ) and 
b = 2nV?T-y 2 : 
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Figs. 3a and 3b depict the dependencies of I* es and 
Rres, on parameters a and b. The resonances can be de- 
tected by peaks in the reflection spectrum, R{y). The 
bright area in Fig. 3a corresponds to well-excited local- 
ized states with high values of the intensity I res , whereas 
the bright area in Fig. 3b indicates easily detectable res- 
onances with pronounced minima in R(v). 
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FIG. 3: The rescaled resonant intensity (a) and reflection 
coefficient (b) via parameters a and 6 

If dissipation is extremely weak so that t « 1, the 
influence of dissipation on the resonances is negligible, 
Ross "C 1, and strongly-excited modes are located near 
the center and in the first, input half of the sample 
(Fig. 3a) as is in a lossless system Q. It follows from 
Fig. 3b that for b <§C 1, easily-detected modes are excited 
only in a relatively small region near the center of the 
sample. Thus, only a relatively small part of excited res- 
onances can be detected via R{v) or T(v) measurements. 

If b > I [but Eq. P still holds], R oss ~ 1 and 
small losses strongly affect the resonance characteristics. 
The absorption substantially reduces resonance transmis- 
sion, T res <C 1, Eq. (3), and also gives rise to the non- 
monotonic intensity distribution inside the first (close to 
the input) effective barrier. The intensity of the reso- 
nances decreases with increasing b (Fig. 3a), and the re- 
gion, where well-excited resonances are located, becomes 
narrower and closer to the input. At the same time, the 
region in which resonant states can be detected becomes 
larger and also shifts towards the input (Fig. 3b). For 
b ^> 1 almost all well-excited resonances are detectable 
via R{v) measurements. Thus, the "detectability" of res- 
onances is improved when small dissipation is present. It 
can be shown from Eqs. IJ7J that the number of reso- 
nances that can be detected in reflection is twice as large 
when 6 > 1 than in the opposite case (b <C 1) for any 
given minimum in the discernible value of R. Moreover, 
at different values of b (or T) resonances in different re- 
gions of the sample are excited and detectable, which can 
provide for the scan of the sample through variations of 
losses. Thus, surprisingly, the presence of loss can make 
it easier to find the characteristics of the eigenmodes of 
a random sample. 

It is important to notice, that for 6^1 the reflection 
coefficient at resonances is a non-monotonic function of 
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dissipation (Fig. 3b), and reaches zero when 
b = a — a -1 . 



(8) 



This effect is well-known in optics and microwave elec- 
tronics as critical coupling [9j • Since no energy is reflected 
from the sample at critical coupling, it corresponds to the 
resonance with the highest intensity I res . Eq. © deter- 
mines the position, x c = L/2 + d, of the most powerful 
resonance in a sample with a given T: 



Hoc 



ln(nVi). 



(9) 



Unlike the lossless case, the disappearance of reflection 
is caused not by the high transparency of the sample 
but rather by strong absorption in the effective cavity. 
The distance well as the maximal intensity does 

not depend on the total length L of the sample. In this 
regard, the sample can be considered as a half-infinite 
random medium. 

Equations ©~© enable one to formulate an algo- 
rithm for retrieval the internal characteristics of local- 
ized states via external measurements. Indeed, quan- 



tities R res , T res , 5v, R, 



l VP 



i — rijocVe — T t 



typ; 



and 



Ttyp — exp (—L/lioc) can, in principle, be measured out- 
side the sample. Then, Eqs. ©-(J!]), constitute a system 
of equations, from which two internal parameters of the 
sample, T and li oc , and three parameters of the reso- 
nances, I, d, and I res can be found. However, even small 
dissipation can make transmittance non-measurable (as 
it was in our experiment), and only reflected signal can 
be analyzed. Nonetheless, the developed model still al- 
lows determination of parameters of the sample and res- 
onances. When b ^> 1 the sample can be treated as a 
half-infinite and one can put Ti = in Eqs. ©-© and 
obtain the following sample parameters: 



r = 



2tt(5i. 



1 - y/%. 
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(10) 



The absorption rate T and the localization length li oc are 
expressed via readily measured quantities R res , 5u and 
Rtyp- The values of T and li oc calculated by means of 
Eq. (|10p with 8u and R res measured for all configurations 
are in agreement within T = 3.76 • 10 _3 cm _1 ± 5% and 
lioc — 1.26cm. The central position, x c , of a resonance 
and its peak intensity, I res , are given by 
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(11) 
(12) 



To check the algorithm of the retreival of the internal 
parameters of resonances from the outside measurements 
of the reflection coefficient, we have compared the coor- 
dinates, x c , and the normalized intensities, I res , of the 
resonances calculated by Eqs. (fTl") (using measured R res 



and 8v) with the corresponding values measured imme- 
diately inside the waveguide by an antenna incerted into 
the waveguide through a slot. As examples, the calcu- 
lated (first line) and measured (second line) values of x c 
and I res for three resonances depicted in Fig. la (^=14.8, 
15.5, 16.5 GHz) are given in the following table: 

x c (cm) 



calc. 4.0, 4.8, 3.6 190, 150, 123 
meas. 4.2, 5.5, 3.6 86, 201, 96 
Note that most of the detected resonances were in 
vicinity of the critical coupling regime (white area in 
Fig. 3b). 

To conclude, localized states in a ID lossy random sys- 
tem has been studied experimentally and analytically. 
The theoretical analysis is based on a model, which de- 
scribes isolated disorder-induced resonances in the frame- 
work of the general theory of resonant systems. It is 
shown that when the dissipation length is much larger 
than the total length of the sample the effect of dissipa- 
tion off resonance is negligible. However, on resonance, 
the effect is dramatic and vanishes only when the dissi- 
pation length is exponentially large as compared to the 
sample length. Otherwise, the energy at a resonant fre- 
quency dissipates substantially due to the exponentially 
high intensity of the field inside effective resonant cavity, 
resulting in the exponential small resonant transmission. 
The reflection coefficient depends non-monotonically on 
the absorption and drops to zero at some values of the 
dissipation length (critical coupling effect). An algorithm 
has been developed for retrieval the decrement of absorp- 
tion and the localization length, as well as the peak in- 
tensities and the positions of the localized modes inside 
a sample from the measurements of the reflection coef- 
ficient. Surprisingly, losses improve the "detectability" 
of resonances and change the conditions for their exci- 
tation. By changing the loss one can scan the sample 
detecting modes excited at different points of the sam- 
ple. The algorithm has been tested by comparison to 
the direct measurements of the microwave field inside a 
disordered single-mode waveguide. 
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